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The measured D*^, D^, and photoproduction cross sections 
have been used to determine charm fragmentation ratios and fractions 
of c quarks hadronising as a particular charm hadron, /(c — * D,Ac). 
Events with a D*^ meson produced in association with an energetic jet 
have been used to measure the charm fragmentation function. The results 
are compared with different models and with previous measurements. 

1 Introduction 

Charm quark production has been extensively studied at HERA using D*± and 
Df mesons. The data have been compared with the theoretical predictions by 
assuming the universality of charm fragmentation. This assumption allows 
the charm fragmentation characteristics, obtained in e^e~ annihilations, to 
be used in calculations of charm production in ep scattering. Measuring the 
charm fragmentation characteristics at HERA permits the verification of the 
charm-fragmentation universality and contributes to the knowledge of charm 
fragmentation. 

The production of D*^, , Df and A^ charm hadrons have been 

measured in ep scattering at HERA in the photoproduction regime {Q^ « 
0) [1] . The measured production cross sections have been used to determine the 
ratio of neutral and charged D meson production rates, Ru/di the strangeness 
suppression factor, 7^, the fraction of D mesons produced in a vector state, 
Pv, and the fractions of c quarks hadronising as a particular charm hadron, 
/(c D,Ac). Events with a D*^ meson produced in association with an 
energetic jet have been used to measure the charm fragmentation function [2]. 
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Figure 1: Distributions of the reconstructed invariant mass for the candi- 
dates (left) and for the candidates (right). The sohd curves represent fits 
to a sum of a modified Gaussian function and a hnear background function. 



2 Measurement of charm fragmentation ratios and fractions 

The production of D*^, D^, D^, Df and charm hadrons was measured in 
the kinematic range pt{D, Ac) > 3.8 GeV and \r]{D,Ac)\ < 1.6. The measure- 
ment was performed for photon-proton centre-of-mass energies in the range 
130 < W < 300 GeV using an integrated luminosity of 79pb~^. Figure 1 
shows distributions of the reconstructed invariant mass for the and Af 
candidates reconstructed from the decay channels K^tt^tt^ (+c.c.) 

and A+ — ^ K^pn^ (-|-c.c.), respectively. The mass distributions were fitted 
to a sum of a "modified" Gaussian function and a linear background function. 
The modified Gaussian function, which was designed for the best description 
of the reconstructed signals, took the form: 



Gauss™°'* (X exp[-0.5 • x^+^+^l 
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Table 1: The fractions of c quarks hadronising as a particular charm hadron. 
The first and second uncertainties represent statistical and systematic uncer- 
tainties, respectively. For the values obtained for charm production in e+e~ 
annihilations, the combined statistical and systematic uncertainties are quoted. 

where x = \[M{Ktt) — Mo]/(t\. The signal position, Mq, and width, a, as well 
as the numbers of charm hadrons in each signal were free parameters of the fit. 
Other details of the charm-hadron reconstruction are discussed in [1]. 
Using the measured cross sections, the charm fragmentation ratios are 

K/d = 1-014 ± 0.068 (stat)t°:°24 (syst), 

7s = 0.266 ± 0.023 (stat)tg;gj^ (syst), 
= 0.554 ± 0.019 (stat)+°;°°^ (syst). 

The measured Ru/d value agrees with one. This confirms isospin invariance 
which suggests u and d quarks are produced equally in charm fragmentation. 
The s quark production is suppressed by a factor « 3.5, as the; measured 

value shows. The measured fraction is sizeably smaller than the naive 
spin counting prediction of 0.75. The predictions of the thermodynamical ap- 
proach [4] and the string fragmentation approach [5], which both predict 2/3 
for the fraction, are closer to, but still above, the measured value. 

The fraction of c quarks hadronising as a particular charm hadron, /(c — > 
D, Ac), is given by the ratio of the production cross section for the hadron to 
the sum of the production cross sections for all charm weakly-decaying ground 
states. The measured fragmentation fractions are compared in Table 1 with the 
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values obtained for charm production in e+e" annihilations [3]. These measure- 
ments, as well as the values obtained in deep inelastic scattering (DIS) [6], agree 
within experimental uncertainties. This confirms the universality of charm frag- 
mentation. 

3 Measurement of chetrm fragmentation function 

Fragmentation fractions are used to parameterise the transfer of the quark's 
energy to a given meson. The measurement of the charm fragmentation func- 
tion in the transition from a charm quark to a D*^ meson was performed for 
photon-proton centre-of-mass energies in the range 130 <W< 280 GeV using 
an integrated luminosity of 120 pb~^. Using events with a D*^ meson pro- 
duced in association with an energetic jet, the fragmentation variable, z, was 
defined as 

where p\\ is the longitudinal momentum of the Z)*^ meson relative to the axis 
of the associated jet of energy E^'^^. The equivalence of {E +p\\)^^^ and 2E^^^ 
arises because the jets are reconstructed as massless objects. The D*^ meson 
was included in the jet-finding procedure and was thereby uniquely associated 
with one jet only. 

Figure 2 shows the normalised differential cross section, 1/ a{da /dz), mea- 
sured in the kinematic range pt{D*^) > 2 GeV, \7]{D*'^)\ < 1.5, E'^* > 
9 GeV and |?7j°'| < 2.4. The data were compared with various fragmentation 
models implemented in the leading-logarithmic Monte Carlo (MC) program 
PYTHIA [7]. The LUND string fragmentation model [8] modified for heavy 
quarks [9] gives a reasonable description of the data [2]. In Fig. 2, the mea- 
surement is compared with PYTHIA predictions obtained using the Peterson 
fragmentation function [10] with different values of the parameter e. The MC 
was fit to the data via a x^-minimisation procedure to determine the best value 
of e. The resuh of the fit is e = 0.064 ± 0.006lo;oJ8- The result is in reasonable 
agreement with the default value used in PYTHIA (0.05), and with the value 
0.053 obtained in the leading-logarithmic fit [11] to the ARGUS data [12]. 
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Figure 2: Relative cross section l/a{da/dz), for the data compared with 
PYTHIA predictions for different values of the parameter e in the Peterson 
fragmentation function. 



4 Summary 

The measured D*^, D", D^, Df and photoproduction cross sections have 
been used to determine charm fragmentation ratios and fractions of c quarks 
hadronising as a particular charm hadron. The measured ratio of neutral and 
charged production rates agrees with one. This confirms isospin invariance 
which suggests u and d quarks are produced equally in charm fragmentation. 
The s quark production is suppressed by a factor w 3.5, as the measured 
value of the strangeness suppression factor shows. The measured fraction of 
D mesons produced in a vector state is sizeably smaller than the naive spin 
counting prediction of 0.75. 

The fragmentation function for D*^ mesons has been measured by requiring 
a jet to be associated with the D*^ meson. The LUND string fragmentation 
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model gives a reasonable description of the data, as does the Peterson function 
with e = 0.064 ± G.006lo;"o8 as determined from a fit to the data. 

All measured fragmentation characteristics agree with those obtained for 
charm production in e+e" annihilations, thus confirming the universality of 
charm fragmentation. 
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